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' The color screening in a QCD plasma, that was studied in a formulation making evident similarities 

' and differences with the electric case, is continued by taking into account the contributions of real 

. gluons. The results, which include a numerical analysis not previously performed, show a damping 

04 ' of the correlation function which, if not exponential, does not differ very much from that form. 

^ ' The role of the temperature, which affect both the population and the dynamics of the quark-gluon 

^ , system, is found to be relevant. 

Q: 

O ! 1. INTRODUCTION 
(N . 

^-H ' The analysis of the screening in gg— plasma performed keeping the analogy with the classical treatment of an electric 

J> \ plasma 0] gave as result a slight difference in the screening behavior with distance, with respect to the more standard, 

r I exponential decay [1,1^^. Now this analysis is pushed further by introducing the effect of the real gluons which 

2^ \ are certainly present in such situations. Once the relevance of the gluons in the quark-quark correlation functions is 

• taken into account one in lead, unavoidably, to consider also the quark-gluon correlations and this is performed in the 
^-H ' present paper. While a definite quark population is taken as a initial condition of the problem, the gluon population is 

' considered of thermal origin and is given therefore in terms of a Bose- Einstein density, together with this population 

\ also the creation of gg— pair is considered, with a corresponding Fermi-Dirac density. Being interested in the large 

1 distance effects, that part of the interaction of real gluons with quark is considered, which correspond to a channel 

exchange of a virtual gluon, so that it has, suitably redefining color and spin factors, the same form of quark-quark 

interaction. 

^ \ The result is a general behavior of the correlation function not very different from that previously found in a 

i-S^ ' pure gg— plasma, since the analytical form is however awkward a numerical investigation has been performed for both 

^ \ situations without gluons and with gluons. The behavior differs little from the exponential, in both case the possibility 
• of small oscillations is found, but only when the correlation function is already very small. 

rS A discussion of the range of values for which the model can be meaningful and also a comparison with some different 

• treatments and results, in particular with the formulation using the termal Green functions are briefly presented j^. 
■ ■ ■ ' Three short Appendices contains some details of the needed calculations. 



2.STUDY OF THE CORRELATION FUNCTIONS 
2.1 General form of the equations 

As it has been explained in the Introduction the basic equation for the correlation functions have the same general 
features as for the pure qq populations, but we must distinguish two possible correlation, the two-fermion case and the 
fermion-boson case. In principle we should look also at the two-boson correlation function, since however our main 
interest is in the two-fermion correlations while keeping into account the boson-fermion correlation, that enter directly 
in the equation for the two-fermion function we neglect the effect of the pure bosonic correlation: in so doing we are 
forced to bound our discussion to situation where the fermion density is dominant, having in mind the nucleus- nucleus 
collision there are certainly situations, of not too large temperature, where this condition is satisfied. 

The ways in which the usual determination of the shielding has been adapted to the case of non commuting charges 
have been presented in the provios paper and will not be discussed here, only the integrodifferential equation for the 
two body correlation is reproduced 
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where the variables q embody both the space variables r and the color indices and the integration over the auxiliary 
variable t takes care of the non commutativity of the charges. 

Now a further derivative with respect to ri is calculated, then one performs the Fourier transform with respect to 
the space variables and the Laplace transform with respect to the inverse temperature with the result: 
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The check (as C) means both the Fourier and the Laplace transform; the variable k is the Fourier-conjugated of r and 
s is the Laplace-conjugated of the inverse temperature (3, Aira = . The factors C{s] k) and T are matrices so their 
order is relevant. In order to proceed is is necessary to specify the colour structure of the quarks and of the gluons, 
so the indices will now must now be displayed and the precise meaning of T and C is specified. It is not convenient 
to work with two set of indices, for quark and gluon respectively, so the color charge of the gluon is represented by a 
traceless tensor with triplet indices, like A'^ with — 0. 

The term T is specified in this way: the interactions qq and qq are given by [App.A] 
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the interaction qq is given by —J; the interaction qg is given by 
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the interaction qg is given by — J. 

The term C is indicated with different symbols following its quark and gluon content and is decomposed according 
to the color structure: for qq and qq it is called respectively Q and Q and decomposed in terms of triplet and sextet 
for qq it is called M and decomposed into singlet and octet and for qg or qg it is called B ov B and decomposed into 
triplet, sextet and 15-plet by means of the following projectors: 
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'Plff = \ 51 5i 5)-\ [Si 51 5) + 51 5l 5j\ + 1 5^ 5^ 5j 
'Paff = \ [5a K si si] Sj 1 [S^ Si S^j. + Si S<i 5j - SI Si Sj Si SI S]] 
''Palcf = I K Si + si si] Sj-^ [Jf si S) + si si Sf + si si Sj + si si S)] (4) 
The correct normalization of the projector is easily verified 
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For the different amplitudes, the equations corresponding to the general form of eq.||2I) are 
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The coefficients p, p, 7 give the densities of quarks, antiquarks ans gluons; the equations for Q and B can be obtained 
by interchange Q Q, p p and B —B in the equations for Q and B. It is possible to give a graphical 
representation to eq.®, see FIG[21 through the rules set in FIG^ it is consistent with the colour structure to use 
lower indices for incoming quarks and outgoing antiquarks and upper indices for incoming antiquarks and outgoing 
quarks. Since the gluon colour structure is given in the spinorial version [App.A], the gluon is represented in the form 
of a qq pair, where up-arrow and down-arrow are used respectively for quarks and antiquarks. 



2.2 Solution of the equations in the conjugate variables 



We recall that the interaction term J is pure imaginary; in the last of eq. © it can be seen that B depends linearly 
on J, in the equations for M and Q, however, B and J appear always together and always multiplied by 7, so it is 
useful to perform the substitutions J = —iJ' ; B = —iB'; in this way in the last equation we get a trivial factor —i, 
in the other two we get a change of sign in the terms that multiplies 7. In the whole set of eq (5) this amounts to 
eliminate the "prime" in J and B and to change the sign 7 — > —7. When we perform the decomposition 
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we find contracting the interaction term J with the projectors P that the RHS of equations ^ for the {qq, qq, qq) 
systems contain only the /^'^ tensor; in fact in any case an octet is exchanged in the t-channel, in this way the relations 
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hold as in the case without gluons and we get the equations: 
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where Bi^ and B15 have been substituted by means of the following relations: 

Bzii^Kit^^^ 3P3+6B6 + 15B15 -0, 3B3+6B6 
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The system ^ yields immediately Fi = F3 + F3 and is reduced to a two-equation system: 
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From the previous linear system it results that F3 and F3 can be expressed in terms of {B^ + Bq) and {B^ + Bq), as 
follows: 
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and n = p + p. Furthermore we get 
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Projecting the equations satisfied by B and B with P the following equations are obtained: 
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The equation for B3, Bq can obtained by the interchange Bi —Bi foi i = 3,6, p p and F3 <^ F3 in the equations 
for B3, Bq. After some tedious calculations, the following two equations for (,63 + Bg) and {B3 + Bq) variable are 
obtained: 
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Substituting the previous expressions in eq.Q we have: 

F - P - , 32^2^2^ 



s [k^ + QkP' SIT an + AStt"^ ^ n + ^'k'^ a? in?) 
32 7r2 a2 ^ 



(10) 



(/c* -\- Q s TT a n + A% t:'^ oP ^ n ~\- S 'k'^ -n?) (fc^ s + 4 tt a n) 
From ea. (|8l9l ?? llU|) it can be easily verified that: 

-83 = ^^3, Bq = ~Bq, Bi5 = —-Bis (11 



2.3 Correlations in real space 



All the gg— correlations may be expressed in term of ^3, so we evaluate its Laplace antitrasform with the following 
result: 
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where y = ■y/487/n — 1 . The inversion of the Fourier transform leads to an expression that is very little transparent, 
it is reported for completeness in [App.C]; the correlation function in real space has been computed numerically, by 
means of standard integration subroutines 0, starting form the following integral 
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In FIG 121 and FIG 0] the correlation function without and with gluons for different choices of temperature, at fixed 
couphng constant a = 0.2 and initial quark density h = 2fm~^, is shown. In both the cases it can be noted that 
the damping of the correlation, if not exponential, differs very little from that behavior; there is also a long distance 
r region (with r ~ 8 /to at 1//3 = 350 MeT^), where the correlation begin to oscillate, but the amplitude of these 
oscillations is very small, being present when the correlation has been already much damped. The previous plots 
give a good overall description of the behavior of the correlation function, but do not show immediately how these 
functions differ from a Yukawa shape e~^^ /r\ in order to give a more complete description of this property in Fig^ 
the expression — ln(rG/3(r)) is plotted. If the Yukawa shape was exact, we would find a straight line; these plots 
confirm (obviously) that the gluons make the damping weaker, moreover the correlation function becomes slightly 
farther from a Yukawa shape. It can be seen that the effect associated to the gluon presence is to produce, at fixed 
temperature, a perceptible increase in the value of the correlation and consequently a displacement towards higher r 
values of the region where the correlation function begin to oscillate. For a detailed analysis of the values of energies, 
quark, anti-quark and gluon densities, we refer to the next section and [App.B]. 



3. CONCLUSIONS 



The inclusion of the gluon in the effect of mutual screening in a gg— plasma refines a previous analysis given in 
term of a pure quark- antiquark population, and it confirms the results. In particular the correlation length and the 
shape of the damping is still the same for qq and qq, so that it must reflect in the same way in the meson and baryon 
production. 

An item that requires some condideration is the comparison with other ways of dealing with the same phenomenon. 
The treatment given in term of strong coupling on the lattice is difficult to compare with |^ |^ because it is very 
different since the beginning, more similar are the treaments in terms of thermal Green function as in In that 
case one extracts the contribution of the pole of the propagator in momentum space so that a precise Yukawa-like 
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decay of the correlation function is certainly produced. The treatment here presented is purely static, so it seems that 
it lacks of some effects that are present in the other treatment, does it contains also someting more? Perhaps what 
is more is seen looking at the structure of the coupled equation ((SJ : the interaction of one particle takes place with 
other particles which are already correlated, so that at least at the level of two-fermion distribution a self consistent 
treatment of the correlations is performed. 

It must be noted that the dependence of the shielding effect on the temperature has two origins: one is the 
kinematical effect that is present in every plasma-like system, the other is dynamical and typical of a relativistic 
system since it comes from the thermal production of particles, both gluons and qq pairs. The two processes have 
not the same role because initially there are more quarks than antiquarks, however, looking at the evolution of the 
chemical potential FIGlHI we see that there is an equivalent temperature, not extremely high, at which the effects 
of the initial condition become very small. One could also ask what the answer would be in the opposite limit, low 
temperature and high density, which might be realized inside a neutron star, but in this situation the Pauli principle 
would be very relevant and one should start from the consideration of a degenerate relativistic plasma, what is not 
attempted in this paper. 

Appendix A: Quark-gluon interaction 

As explained previously only the quark-gluon interaction corresponding to a gluon exchange in channel is considered. 
For this amplitude tha color structure is usually given in mixed form [|| as fABc{Tc)1 where A,B,C = 1 ... 8 and 
a, 6 = 1,2,3. This form is not convenient here and will be transformed in a purely spinorial version. Using the 
fundamental commutator [Ta,Tb] = ifABcTc we express the interaction as commutator; then defining the matrices 
with spinorial indices 

the form of the interactions used in this paper for qq and for qg is 
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Note that the antisymmetry of J with respect to the exchange (c, d) ^ (/, g) together with the factor i ensures the 
Hermiticity of the interaction. 



Appendix B: Densities of quarks, antiquarks and gluons 

When we take care of the fact that there can be a thermal production of gluons we cannot ignore the concurrent 
production of quark-antiquark pairs. We investigate briefly the problem in this form. Assume the baryonic density, 
i.e. the quark density minus the antiquark density, as given and find the simultaneous production of gluons and of 
quark-antiquark pairs. The expression for the baryonic density is simple when the quark mass can be neglected, or 
gives at most a small correction, in this case for the density we get: 
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In particular the antiquark density is given by 
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where Ck{z) = X^^i ^"Z*^*^ ^^^d for the weight we get gf — 12 resulting from a factor 2 for the spin, a factor 3 for 
the color, a factor 2 for the flavor because we consider only u and d quarks. 
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From eq (|Tl|l it follows: 
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From this relation we see the conditions qualitatively well known in which the mass term is negligible: either m << l/j3 
or m << /i, form now on we shall assume that at least one of these conditions is fulfilled and the mass term is dropped. 
Solving the equation for we find 
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with X = bf3^. In this way we get the actual expression of p in terms of b and /3, the total density which appears in 
the previous formulae is given by n = 5 + 2p. 

It is now possible to give a semi-quantitative analysis of the different densities which are relevant for our problem: 
we start assigning a density b which we take equal to 2 fm^"^. In FIGEI FIGEjand FIG|Hlthe values of /3fi, n and the 
ratio "f/n as functions of 1/(3 are shown. 

We have worked out three cases: the first is chosen so that the density of antiquarks is much less than 6, we take 
1/(3 = 150 Mel/, this gives P ~ -^b and for the corresponding gluon density 7/n = 0.40, in the second case we take 
1/(3 — 300 MeV, this gives p « | 5 and 7/n = 0.84, in the third one 1/(3 — 450 MeV is considered, which produces 
p K, Qb and 7/n — 0.88, which is the "saturation" value for 7/n ratio. In these expression the density of gluons has 
be taken as the free-boson thermal density, which amounts to: 



(18) 



Here gi, = 16 i.e. a factor 2 for the spin and a factor 8 for the color. Since we have (3p — > for /? ^ it results that 
in that limit 
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and this limit is numerically reached already at 1/(3 — 390 MeV. 

Appendix C: Analytical expression of the correlations 

The inversion of the Fourier transform 
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where the generalized hypergeometric functions are used 
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and the coefficient are 
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The function Gp{r'^) is real because it is the sum of two terms with their complex conjugate. Note that for very low 
values of 1/(3, y could become imaginary, in this case all the addenda would be separately real. 
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FIG. 4: Correlation function Gg{r) in presence of gluons for different choices of temperatures, at fixed coupling constant a = 0.2 
and initial quark density b = 2fm~^. 
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a = 0.2 b = 2 fm"^ 1/p = 150 MeV 

: with gluons ^ ' 

■ : no gluons 



a = 0.2 b = 2fm"' l/p = 300MeV 

: with gluons 

■ : no gluons 




a = 0.2 b = 2fm"^ l/p = 450MeV 

: with gluons 

■ : no gluons 



FIG. 5: Minus Logaxithm of rG/3(r) with (solid line) and without (dashed line) curve for different choices of temperature, 
coupling constant a = 0.2 and initial quark density 6 = 2 fm~^. 
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FIG. 6: P n as & function of with initial quark density 6 = 2 fm ^; note that it is an adimensional quantity. 
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FIG. 8: 7/n ratio as a function of 1/(5 with initial quark density given by 6 = 2 frri 



